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Abstract. The solution Xn to a nonlinear stochastic differential equa- 

1 Y^jv 

2 ■^j = 



tion of the form dX,,{t) + A,^{t)X„{t) dt ~ i J2'^^iiBjit)fXn{t)dt 



EjLi B;i(t)X„{t)dpj^(t) + fn{t)dt, Xn{0) = X, where is a regular ap- 
proximation of a Brownian motion /3j, B"{t) is a family of linear contin- 
uous operators from V to H strongly convergent to Bj{t), A„{t) A{t), 
{Ari{t)} is a family of maximal monotone nonlinear operators of subgra- 
dient type from V to V' , is convergent to the solution to the stochas- 
tic differential equation dX{t) + A{t)X{t)dt - \ Y!^^^ B]{t)X{t) dt = 
Y.J=iBj{t)X{t)dl3j{t) + fit)dt, X{0) = X. Here V C H ^ H' C V 
where is a reflexive Banach space with dual V' and H is a Hilbert 
space. These results can be reformulated in terms of Stratonovich sto- 
chastic equation dY (t) + A{t)Y (t) dt = Bj{t)Y{t)odPj{t) + f{t)dt. 
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1. Introduction 
Consider the stochastic differential equation 

( N 

dX{t) + A{t)X{t) dt-lY^ B]{t)X{t) dt = 

TV 

Bj{t)X{t) d(3j{t) + fit) dt, t € [0, T] 

i=i 

I X{0) = X 

(1.1) 

where A{t): V ^ V' is a, nonlinear monotone operator, Bj[t) € L{V,H), 
G [0, r] and f5j are independent Brownian motion in a probability space 

{n,^,{^t}t>o,n- 

Equation (jl.ip is of course equivalent with the Stratonovich stochastic 
differential equation 

(1.1)' dY{t) + A{t)Y{t) dt = ^f^-^Bj{t)Y{t) o df3j{t) + fit) dt, te [0,T]. 
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Here ^ is a reflexive Banach space with dual V such that V C H C V 
algebraically and topologically, where H (the pivot space) is a real Hilbert 
space. (The assumptions on A{t), Bj{t) will made precise later on.) 
We associate with (jl.ip the random differential equation 

-£+A{t)y = git), t E [0,r], P-a.s. 

(1.2) 

. y(0) = X, 

where g{t) = e^^''=i ''•'^*^^^*(*^/(t), B*{t) is the adjoint of Bj{t) and A(t) is 
the family of operators 

A(i)y = e-Ef.iftWB,W^(^)gEf.,ft(t)B.(t)y + 

+ V / ' e-^-^^W^j(t)e^-^^Wyds, Vy G F (1.3) 

where Bj is the derivative of t — ^ Bj{t) G L{V,H) and (e^-^j(*^)seR is the 
Co-group generated by Bj{t) on H and y. 

It is well known (see e.g., [T71 pag. 202], [H], [21] ) that assuming that 
BjBk{t) = Bk{t)Bj{t) for all j, fc, at least formally equations (jl.ip and (|1.2p 
are equivalent via the transformation 

X{t) = eSf=i^^W^^(*)y(t), P-a.s, t € [0,r], (1.4) 

and this is indeed the case if p.2p has a strong, progressively measurable 
solution y : [0, T] X ^ if. 

We consider also the family of approximating stochastic equations 

d ^ 

-Xn + An{t)Xn = B^ {t)Xn{t) $J (t) + /„(t), P-a.S. 

j=i (1.5) 

XniO) = X. 

where {Pj"} is a sequence of smooth stochastic processes convergent to (3j, 
that is P^{t) I3j{t) uniformly on [0,r], P-a.s. and An A, Bf Bj, 
/n — > / as n — > oo in a sense to be made precise below. 

Equation (jl.Sp is just an approximation of Stratonovich equation (jl.ip ^ 
where is a regularization of /3j. One must emphasize that {/3n} might be 
adapted to a filtration different of {^t}- 
Equation (jl.Sp reduces via (jl.4p . that is 

(1.4)„, X„(t) = eS.^^/^"W^"(Sn(t) 

to a random differential equation of the form (|1.2p that is 



(1.2)„, 



^ + K{t)yn = gn{t), t G [o,r], p-a.s. 

yn(0) = X 
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The main result (see Theorems 12.11 12.21 12.31 below) is that under suitable 
assumptions equations ()1.2p . (1.2)„ have unique solutions y and y„ which 
are progressively measurable processes and for n — )■ oo, we have yn — > y, 
Xn — > X in a certain precise sense. In the linear case such a result was 
established by a different method in [T?] , (we refer also to [18] , [21] to other 
results in this direction.) The variational method we use here allows to treat 
a general class of nonlinear equations (jl.ip possibly multi-valued (On these 
lines see also [S1[B].) 

The applications given in Sect. 4.1 refer to stochastic porous media equa- 
tions and nonlinear stochastic diffusion equations of divergence type but of 
course the potential area of applications is much larger. 
Notation. If y is a Banach space we denote by L^(0, T; Y), 1 < p < oo the 
space of all (equivalence classes of) F-measurable functions u: (0,T) Y 
with \\u\\y e LP(0,r) (here || • ||y is the norm oiY). Denote by C{[0,T];Y) 
the space of all continuous valued functions on [0, T] and by T4^^'^([0, T]; Y) 
the infinite dimensional Sobolev space {y € L^{0,T;Y), ^ G LP{0,T;Y)} 
where ^ is considered in sense of vectorial distributions. It is well known 
that W^'P{[0,T];Y) coincides with the space of absolutely continuous func- 
tions y: [0,r] — > Y, a.e. differentiable and with derivative y'{t) = ^{t) a.e. 
t e (0,T) and ^ G LP{0,T;Y){see e.g., [1].) If p e [l,oo] is given we denote 
by p' the conjugate exponent, i.e., ^ + ^ = 1- 



We shall study here equation ()1.2p under the following assumptions 

(i) V is a separable real reflexive Banach space with the dual V and H 
is a separable real Hilbert space such that V C H C V algebraically and 
topologically. 

We denote by | • |, || ■ ||v^ and || • ||y the norms in H, V and V, respectively 
and by (•, •) the duality pairing on V x V which coincides with the scalar 
product (•, ■) of H on H X H. 

(ii) A{t)y = d%l){t,y), a.e. t G (0,r), Vy G V, F-a.s., where ip: {0,T)xV x 
Q M is convex and lower- semicontinuous in y on V and measurable in t 
on [0, T]. There are > 0, 7^ G M, i = 1, 2, 1 < pi < p2 < og 



2. The main results 



71 + ai \\y\\v 



< i^{t,y) < 72 + "2 \\y 



\P2 

\v ' 



Vy G V. 



(2.1) 



(iii) There are Ci, C2 G such that 

m-y) <Ci^{t,y) + C2, VyGKtG (0,r). (2.2) 
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(The constants Cj, 7^, Ui are dependent on a;.) 

(iv) For each y G V the stochastic process ipi^t, y) is progressively measurable 
with respect to filtration {^t}t>o- 

( v) Bj [t) is a family of linear, closed and densely defined operators in H such 
that Bj{t) = -B*{t), Vi e [0,T], Bj{t) generates a Co-group (e'^-^^Wj^gR on 
H and V. Moreover, Bj G C^{[{),T]-L{V,H)), Bj{t)Bk{t) = Bk{t)Bj{t) for 
all j, k. 

(vi) /: [0, T] X — )• y' is progressively measurable and f € Lp'i{0,T;V'), 
¥-a.s. 

We note that by (ii) A{t,u}): V ^ V is, for all t G [0,T] and w G 17, 
maximal monotone and surjective (see [3]) but in general multi- valued if ip 
is not Gateaux differentiable in y. 

Theorem 2.1. Let y^ ^ H . Then under assumptions (i) ~ (vi) there is for 
each cj G a unique function y = y{t,uj) to equation (|1.2p which satisfies 

y G LPi(0,r;y) nC([0,r];i7) n VF^'P2([0,r];T/')> (2-3) 

'^{t)+K{t)y{t)^ g{t), a.e.t£{0,T), 
'^^ (2.4) 

y(0) = yo- 

Moreover, the process y : [0, T] x J7 — >■ if is progressively measurable with 
respect to the filtration {^t}t>o- 

If An{t, •) and A^^{t, •) are multi-valued we mean by An{t, y) and A~^{t, z) 

single valued sections. 

G 

By — > we denote the variational or F-convergence. This means that for each 
y & V and ^ G A{t)y there are y„ and ^„ G A{t)y such that Vn ^ y strongly 
in V , £,n ^ ^ strongly in V' and similarly for A~^{t) — > A~'^{t). Assumption 
(j2.5p implies and is equivalent to: ipnit,z) — > ip{t,z), tp^{t,z) — > ip*{t,z) for 
all z & V, z £ V and t G [0, T] where V'* is the conjugate of ip (See e.g.. 



Theorem 2.2. Assume that for each n, A„, i?" and /„ satisfy (i) ~ (iv). 
T/ien /or any yg G F i/iere is a unique function yn = yn{t,oj) which satisfies 
(|2.3p and equation (|2.4p iwii/i A„ instead of A.. Moreover, assume that for 
n — )• 00 

^„(t) 4 A(t), t G [o,r] 
A-i(t)^^-Ht), tG[o,r] 

/„(-,a;) ^ /(-.w), strongly in LP'^{Q,T-V'),¥-a.s. inn. (2.6) 
B^x^BjX, inC\[0,T];H), e H. (2.7) 

T/ien /or n — )• 00 

y„(-,a;) ^ y(-,a;) (2.8) 
P-a.s. weakly in L^^ (0, T; V), weakly-star in L°°(0, T; if). 
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Assumption (j2.5p implies and is equivalent to: ipn{t, z) — > z)^ V'n(^) z) — > 
ip*{t,z) for all z ^V, z e V and t G [0,7"] where V'* is the conjugate of 
(See e.g., 0.) 

Coming back to equation (jl.ip we say that the process X : [0, T] — )• 
is a solution to (jl.ip . if it is progressively measurable with respect to the 
filtration {^t}t>o induced by the Brownian motion, 

X eC{[0,T],H)r\LP'{0,T;V), P-a.s. (2.9) 

and 

rt 1 ^ 

X{t) = X- A{s)X{s)ds + -y" / B^{s)X{s)ds (2.10) 

+ V / Bj{s)X[s)dpj{s)+ j f{s)ds, VtG[0,r], P-a.s.. 
j^i Jo Jo 

By Theorem 12.11 and Theorem 12.21 we find that 

Theorem 2.3. Under the assumptions of Theorem \2.1\ there exist unique 
solutions X and Xn to (jl.ip and (|1.5p respectively given by 

X(t) = e^." ^^^■W^^(S(t), X^t) = e^^^^^"W^"(Sn(t), (2.11) 

where y and yn are solutions to ()1.2p and (|1.2p ^ . Moreover, we have 

X,Xne LP^{0,T;V), F-a.s. (2.12) 

X,Xn: [0,T] —7- H are P-a.s. continuous and 

Xn^ X weakly in {0,T;V), weakly-star in L°°{0,T;H), P-a.s. 

(2.13) 

The precise meaning of Theorem l2.2l and Theorem l2.3l is the structural sta- 
bility of the Ito stochastic differential equation (jl.ip and of its Stratonovich 
counterpart (1.1)'- As a matter of fact, as mentioned earlier, all these results 
can be reformulated in terms of Stratonovich equation (1.1)'. 
One of the main consequences of Theorem 12.21 is that the Stratonovich sto- 
chastic equation is stable with respect to smooth approximations of the 
process B{t)Xd/3{t). On the other hand, the general existence theory for in- 
finite dimensional stochastic differential equations with linear multiplicative 
noise (see e.g., [I6l[l7]) is not applicable in present situation due to the fact 
that the noise coefficient x — )> B{t)x is not bounded on the basic space H. 
The approach we use here to treat equation (j2.4p relies on Brezis-Ekeland 
variational principle [12], [13j which allows to reduce nonlinear evolution 
equations of potential type to convex optimization problems. (On these 
lines see also l5]-[7], [22], [M])- 

The more general case of nonlinear monotone and demicontinuous operators 
A{t) : V ^ V' is ruled out from present theory and might expect however 
to extend the theory to this general case by using Fitzpatrick function for- 
malism (see [22], [21]). 

As in [Uj, see Corollary on p. 438, we can define a solution to problem 
(1.1)' for any deterministic continuous function /3 : M+ M^. The result 
from jl4] was a generalisation of a analogous result from Sussmann's well 
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known paper [23], see also Doss [19]. We will formulate our result in the 
same fashion as in [TJ], i.e. the result contains implicitly a definition. Let us 
observe, that in this case, we prove the existence for any multidimensional 
continuous signal, thus a signal more general than a rough signal from the 
theory of rough paths. However, this is due to the assumption of the com- 
mutativity of the vector fields Bj, j = I, - ■ ■ ,N. In the result below, we need 
deterministic versions of assumptions (ii), (vi). Note that the assumption 
(iv) is now redundant. 

(ii') A{t)y = di}{t,y), a.e. t G (0,r), Vy E V, where -0 : (0, T) x F ^ M 
is convex and lower-semicontinuous iny onV and measurable in t on [0, T] . 
There exist ai > 0, 7^ € M, i = 1, 2, I < pi < P2 < 00, such that 

71 + ai l|y lly <V'(i,y) < 72 + 02 lll/ll^', G (2.14) 
(vi) f eLPiiO,T-V'), 

Theorem 2.4. Assume that the assumptions (i), (iii), (v) as well as (ii') 
and (vi') are satisfied. Then for every x & V and every j3 G C([0,T];R^), 
the problem 

iN 
dX{t)+A{t)X{t)dt = ^ Bj{t)X{t) df3j{t) + f{t) dt, t G [0,T] 
X{0) = X 

(2.15) 

has a unique solution X G L^^i^O, T; V) PI C([0, T]; H) in the following sense. 

(i) For every (3 G C-'^([0, T]; R^), the problem (|2.15|) has a unique solution 
X G LP^{0,T;V)nC{[0,T];H). 

(ii) If(3n G C\[0,T];R^) and (3n ^ /? in C{[0,T];R^) and Xn G LP^iO,T;V)n 
C{[0,T]; H) is the (unique) solution to the problem ()2.15p corresponding to 
Pn, then Xn X in weakly in L^^ (0, T; V), weakly-star in L°°(0, T; H), P-a.s. 

From Theorem 12.41 we infer that in the framework of Theorem 12.41 but 
with /? being a Brownian motion, the problem (jl.ip generates a random 
dynamical system on H. In an obvious way we have the following Corollary 

Corollary 1. Assume that the assumptions (i), (iii), (v) as well as (ii') 
and (vi') are satisfied. Assume that P in a standard canonical two- sided - 
valued Brownian motion on a filtered probability space {fi, {^(}t>o, P}, 
where = {w G C(M,R^) : a;(0) = 0}. Let us define a map 

{} -.mxnB {t,uj)^i^toj = u}{- + t)- w(o) G n. 

Then there exists a map 

If : X Q X H 3 {t,ijj,x) t-^ ip{t,ijj)x G H 

such that a pair ((j), •&) is a random dynamical system on H , see for instance 
Definition 2.1 in [K], and, for each s G M and each x € H, the process X, 
defined for a; G and t > s as 



X{t, s; uj, x) := ip{t — s; 'i}s^)x, 



(2.16) 
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is a solution to problem (jl.ip over the time interval [s, oo) with an initial 
data given at time s. 

Remark 2.1. Theorem 12.41 provides a solution to problem ()2.15|) for every 
continuous path. Our main result provides a natural interpretation of this 
solution in the case when /3 is a Brownian motion. One can also provide a 
similar interpretation when /3 is a fractional, see for instance |20j . 
Corollary allows one to investigate the existence of random attractors, see 

m- 

These questions will be investigated in the future works. 

3. PROOFS 

Proof of Theorem 2.1 
For simplicity we consider the case = 1, that is Bj = B, f3j = (3 for all j. 
We note first that though the operator A(i) = A(t,a;) is P-a.e w € il, 
maximal monotone from V to V' the standard existence theory (see e.g., 
[H pag. 177]) does not apply here. This is, however, due to the general 
growth condition ()2.ip on '(/'(i, •) and implicitly on A[t) as well as due to the 
multivaluedness of A{t). 

So we shall use a direct approach which makes use of the variational structure 
of problem (|1.2|) . (On these lines see also [S], [71 page 280]). Namely, we can 
write 

K{t) = d^{t,-)+T{t), VtG[0,r]. (3.1) 
Here : [0, T] x F ^> M is given by 

(^(t,2/)=V(t,e'^(*)^«y) (3.2) 

and 

T{t)y= e-^^W^(i)e^^Wyds, Vy € t € [0, T] 

JO 

where B = %B{t). We fix w G 

By the conjugacy formulae (|A.3p and ()A.4p we now may equivalently write 
(fOjl (or dM])) as 

/ ^(t, y{t)) + ip*{t, u{t)) = {y{t),u{t)), a.e. t G [0, T] 

\ y'{t) + T{t)yit) = -u{t) + g{t), a.e. t G [0, T] ^''■''> 

while 

'/'(*, y) + 'P*{t,u) > {y,u) 

for all G LPi(0,T;F) x LP2(0,r;y')- 

Thus following a well known idea due to Brezis and Ekeland (see e.g., 
[?. 112^ [T3]) we are lead to the optimization problem 

Min{ f {^{t, yit)) + ^*{t, u{t)) - {u{t),y{t))) dt : 
Jo 

y + r(t)y = -« + £/, a.e. t G [0, T]; y(0) = yo, 

y G (0, T; y), n G (0, T; )}. (3.4) 
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Equivalently 

rT 



Min{ / {ip{t,y{t)) + ^*{t,u{t)) - {g{t),y{t))) + i(|y(r)|2 - jyoP) : 
Jo 

y' + T{t)y = -u + g, a.e. t € [0, T]; y{0) = yo, 

y G LP' (0, T; V),u G Lp'-{0, T; V')}. (3.5) 

Here we have used (for the moment, formally) the integration by parts for- 
mula 

"^{u{t),y{t))dt = l{\y{T)\'-\yo\') 



{Tit)y{t),yit))dt- / {g{t),u{t))dt 



T 





and hypothesis (v) which implies that {T{t)y,y) = 0. Of course the equiva- 
lence of ()3.4p and ()3.5p is valid only if the above equality is true which is not 
always the case in absence of some additional properties of minimizer y to 
allow integration by parts in {u{t), g{t)) dt. In the following we shall prove 
however that equation ()3.5p has at least one solution and show consequently 
that it is also a solution to equation (j2.4p . 

Lemma 3.1. There is a solution y* G LP'{0,T;V) n W'^'P'2{[0,T];V') to 
equation (|3.5p . 

Proof. We note that by standard existence theory of linear evolution equa- 
tions for each u G LP2 (0, T; V') there is a unique solution y G LP'^ (0, T; V) fl 
W^'P'^{[0,T];V') C C{[0,T];H) to equation 

y' + r{t,y) = -u + g, a.e. tG[0,T], y{0) = yo- 

By assumptions (|2.ip and (|2.2p we have 

7i+ai||y||y <¥'(i,y) <72 + a2||y||^', Vy G F (3.6) 

and 

li+ai\\yfy^<^*{t,y)<j2 + ai\\y\fy\, Vy G (3.7) 

where 7,, 7^ G M, ^ + A = 1 and Q!j,Q;i > 0, i = 1,2. (Recall that e'*'^^*) is 
invertible) . 

Then the infimum m* in (|3.5p is > —00 and there are the sequences {yj} C 
LPi(0,r; V), {uj} C LP2(0,r; F') such that for ah y 

m*< / (((/p(t,y,-) + ^j) - rfi + i(|yi(r)l' - |yo|') < "I* + - 

Jo J 

(3.8) 

y' + T{t)yj = -Uj + g, in [0, T] 

vm = yo. ^'-'^ 

Clearly yj G VF1'P2([0, T]; F') and by assumption ([XT]) and inequality ([^ 
it follows that 

l|yillLPi(o,T;y) + llyjlliP?,(o,T;y') ^ ^ (^-lO) 



because as easily seen by assumption (v), \T(t)y\H < CHyHy, Vy G V. 
Hence on a subsequence, again denoted yj, we have for j — )• oo 

Vj y* weakly in L^i (0, T; V) 

Uj u* weakly in Lp'2{0,T; V) (3.11) 
y'. {y*y = g-u* - T{t)y* Weakly in ^(0, T; V). 

Since the functions y i— )■ ip{t,y) dt and u i— )• Jq ip* {t,u) dt are weakly 
lower-semicontinuous on L^'i(0,r;y) and LP2(0,T;y) respectively, letting 
j tend to infinity we obtain that 

m* = r {^it,y*) + ^*it,u*) - {g,y*)) dt + l{\y*iT)\' - |yoP) (3.12) 



and 

{y*y + rit)y* = -u* + g, t e [0,T] 

y*{o) = yo. ^'^■^^> 

Therefore {y* ,u*) is a solution to optimization problem ()3.5p as claimed. □ 

Proof of Theorem 2.1 ( continued). We shall show now that y* given by Lemma 
I3.1l is a solution to ()2.4p . To this end we notice just that without any loss of 
generality we may assume that yo = 0. Indeed we can reduce the problem 
to this case by translating in problem 12.31 ?/ in y — yo- 

We prove now that m* = 0. For this purpose we invoke a standard duality 
result for infinite dimensional convex optimal control problems, essentially 
due to R.T. Rockafeller. Namely, one has (see [TJ Thm. 4.6, pag. 287]) 

m* + min ([33])' = (3.14) 

where (|3.5p ' is the dual control problem 

Mm{j^ ((^(t, -p{t)) + ^\t,v{t)) + {9{t),p{t))) dt + i|p(r)|2 : 

p' + mp = v + 9, te[o,T]}. 

If ip*,v*) G LPi(0,T;y) X Lp'2{0,T;V') is optimal in (3.5)', we have 

{{pn',p*)^L\0,T) (3.15) 

((p*)',p*)di = i(b*(r)|2-|/(0)|2). (3.16) 

Here is the argument. First, note that p' solves p' + T{t)p = v + g. We have 
by the identities (jX3]) and (jXi]) and the fact that {T{t)p*,p*) = 0, 

-{{p*{t))',p*{t))<^*{ty{t))+^{t,-p*{t))-{g{t),p*{t)), a.e. tG(0,r) 

and 

{{p*{t))\p*{t))<^*{t,v*{t))+^{t,p*{t)) + {g{t),p*{t)), a.e. tG (0,r). 

Since ip{t,-p*) G L'^{0,T) and by assumption ([22]), v?(t,p*) G Li(0,T) too, 
we infer that (|3T5]l holds. Now since p* G W'^'P'2{[0,T];V') D LP^{0,T;V') 
we have 

^^|P*WP = ((/)'W,/(0), a.e. tG(0,r) 
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and by (j3.15p we get (j3.16p as claimed. 
By (3.5)' and (|3.16p we see that 

[■T 

min(3.5)'= / {^{t, -p*) + ^* {t,v*) + {v\p*)) dt + \\p* {fd)\'^ >{). 
Jo 

Similarly, by (IHTTD and by 

i(|y*(r)p-|y*(0)p)= r {{y*)\y*) dt, 

Jo 

(the latter follows exactly as ()3.16p ) we see that 

T 

{ip{t,y*) + ^*{t,u*)-{u\y*))dt>0. 

Then by ()3.14p we have that m* = and therefore again by ()3.12p we have 
that 

ip{t,y*) + ip*{t,u*) = {u*,y*), a.e. in [0, T] 
{y*y + r{t)y*=g-u* a.e. in [0, T] 

and therefore y* is a solution to ()2.4p . 

On the other hand, as seen earlier, we have 

l{\y*{t)\^-\y*is)\^)= [\{y*nr),y*{T))dr, \fO<s<t<T. (3.17) 

J s 

Hence y* € C{[0,T]; H). The uniqueness of y* is immediate by (j3.17p . It 
remains to be proven that y* is progressively measurable. 
To this end we note as minimum in (j3.5p the pair {y*,u*) is the solution to 
Euler-Lagrange system (see e.g., [3 page 263]) 

{y*y + T{t)y* = -u* + g, a.e. t G (0, T), w G Q 
q' - T'{t)q = -g + A{t)y*, a.e. t G (0, T) 
u*{t) = A{t)q{t), a.e. t G (0,r), w G O 
2/*(0)=yo, q{T) = -y*{T). 

Since the latter two point boundary value problem has a unique solution 
{y*,q) and is of dissipative (accretive) type it can be solved by iteration or 
more precisely by a gradient algorithm (see [TJ page 252]). In particular, 
we have y* = liuik^^yk, q = limfc_^oo 9fc weakly in LP^{0,T;V) and u* = 
lim^._^oo Uk weakly in LP'^{0,T; V) where 

yi + r{t)yk = -uk + g t e [0,T], 
q'^-T'{t)qk = -g + A{t)yk, t e [0,T], 
Uk+i = Uk-A ^{t)uk + qk, t£[0,T], 
yk{0)=yo, qk{T)=0, fc = 0,1,2,.... 

Hence, if we start with a progressively measurable uq, we see that all Uk are 
progressively measurable and so are u* and y*. □ 




Proof of Theorem 2.2. As in the previous case it follows that equation ()1.2p „ 
has a unique solution y„ G VF-'^'P2([0, T]; V') PI ^^^(0, T; V) given by the min- 
imization problem (j3.5p where g = g-n and = (^* = Here is 
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given as in (|3.2p where ip = ipn and /? is replaced by while 99* is the 
conjugate of (fn- We have, similarly, dipn = and (pn{t,y) = ipn{t,e^"^^^y) 



dy„ 



dt 



L^2 {0,T;V') 



+ \yn{T)\<C, 



{Vn, Un) = arg min{ / ipn{t, y{t)) + ip*^{t, u{t)) - {gn{t), u{t)) dt 
Jo 

+ i(|y(r)|2-|yo|'); y' + Tny = -u + gn, y{0) = yo}- 

Here r„(t)y = /g^"^*^ e-^-^"WS„,(i)e^-^"Wy ds. 
We see that 

\\yn\\L^{0,T;H) + lll/n || LPi (0,T;y) + 

and this implies that on a subsequence, again denoted {n}, we have 
' Un^u weakly in LP2 (0, T; V) 

yn — > y weakly in L^i (0, T; V) 

yn — > y weakly-star in L°°{0,T; H) 

dyn dy 



(3.19) 



dt 



dt 



weakly in LP^{0,T;V') 



. yn{T) — > y{T) weakly in H. 

By ([23]), (fTB]) we see that y' + Ty = -u + g. 
Moreover, we have 

{iPn{t,yn{t))+ip*^{t,Un{t)) - (5n(t), ?/n(t))) + i |y„,(r) p 



< / {Mt,y*{t)) + ^l{t,u*{t))-{gr,{t),u*{t)))dt+l\y*{T)\^. 
Jo 

where {y*,u*) is the solution to (jS.Sp . Now by assumptions ()2.5p and ()2.6p 
we have 

^*^it,u*„{t))^^it,u*{t)) vtG[o,r], 

and this yields 

limsup / {(pn{t, yn) + Vnit, Un) - {gn,yn))dt + l\y^{T)\^ -l\yo\^ = 0. 

n— >oo JO 

(3.20) 

In order to pass to limit in ()3.20p we shall use ()3.19p and the convergence 
of {(fn} and {99*} mentioned above. We set z{t) = e''^*) ^(*)y(t), z„(t) = 



B{t) 



yn{t)- We have 



Mt, S{t)) < Mt, Zn{t)) + {An{t, S{t)), Z{t) - Zn{t)), 

and since dip^ = we have also that 
rnit, 0{t)) < rnit, On{t)) + {A-\t){t, 9{t)),9{t) - 9n{t)) , a.e. t G (0, T). 



WOREL BARBU, ZDZISLAW BRZEZNIAK, ERIKA HAUSENBLAS, AND LUCIANO TUBARO 

where e{t) = g{t) - y'{t), Onit) = gn{t) - y'^{t). 

Then by assumption (j2.5p and equation (j3.20p we have that 

hmsup / {^n{t,y{t)) + Vnit.9{t)-y'{t))- 

- {g{t),y'{t)))dt+'^\r{T)f -k\yo\' = 0, 

and so, since as seen eadier (j2.5p imphes that ipn{t,z) — >■ Lp{t,z), \/z G V, 
(fnit, z*) ip*{t, z*), Mz* G V\ by the Fatou lemma, we have 

r {ip{t,y)+^*{t,u) - {g,y)) dt + \\y*{T)\^ - \\yo? < 0, 
Jo 

which imphes as in the previous case that y is a solution to ()3.5p and therefore 
to ([2^ as claimed. □ 

4. Examples 

The specific examples to be presented below refer to nonlinear parabolic 
stochastic equations which can be written in the abstract form (jl.ip where 
A{t) are subpotential monotone and continuous operators from a separable 
Banach space V to its dual V'. 

We briefly present below a few stochastic PDE to which the above theorems 
apply. We use here the standard notations for spaces of integrable functions 
and Sobolev spaces wd'^(^), l^-^'^V) = (W^o '^(^))', ^^(^), A; = 1,2 on 
open domains C W^. 

4.1. Nonlinear stochastic diffusion equations. Consider the stochastic 
equation 

' dXt - divg a{t, V^Xt) dt - ^ b{t, C) • V^{b{t, C) ■ V^Xj) dt 

= b{t,0-V^Xtd(3{t), in(0,T)x^ 
Xo = xmff ^ ' ' 

^ Xt = on (0, T) X dff 

Here a: (0, T) x M'^ — > R'^ is a map of gradient type, i.e., 

a{t,y)=dj{t,y), VyGM^tG(0,r) 

where j : (0, T) x M"^ x ^ M is convex in y, progressively measurable in 
{t, w) G [0, T)xn and 

7i + ai|yr <j(i,y) <72 + a2|y^^ Vy G w G fi, i G (0, T) (4.2) 

j{t,-y) <cij{t,y)+C2, Vy GM^t G (0,T). (4.3) 
It should be emphasized that the mapping r — )■ a{t, r) might be multivalued 
and discontinuous. As a matter of fact if o(t, •) is discontinuous at r = rj, 
but left and right continuous (as happens by monotonicity) it is replaced by 
a multivalued maximal monotone mapping d obtained by filling the jumps 
at r = rj. 

Equation (jH]) is of the form 1^ where H = L'^{^), V = Wq'^'{^), 
A{t) = d^{t, •), 2 < pi < p2 < oo, 

i^{t,u) = I j{t,Vu)d^, W G <'^^(^) 
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and 

B{t)u = b{t, • V^u = div^{b{t, C)u), Vu G W^'P' (^). (4.4) 
As regards the function b{t, r) : [0, T] x M'^ — )• we assume that 

h{t,-l |^(i,.)G(C([0,T];^))'^ (4.5) 

r b{t, ■) + ar is monotone for some a > 0, (4-6) 



div^6(i, 0) = 0, b{t, • u{C) = Ve G (4.7) 

where ly is the normal to dff. (The boundary dff is assumed to be of class 
CK) 

Here div^6 is taken in the sense of distributions on 0'. 

Then (4.4) defines a linear continuous operator B{t) from V to H = L'^(ff) 
which as early seen is densely defined skew-symmetric, that is —B{t) C B*{t) 
Vt G [0,T]. Moreover, B{t) is m-dissipative in L^(^), that is the range of 
u ^ u — B{t)u is all of L^(^). Indeed for each / G L^(^) the equation 
u — B[t)u = f has the solution 

/•oo 

"(0= / e~'f{Z{s,C))ds, V^G^, 
Jo 

where s Z{s, is the differential flow defined by equation 
dZ 

— = b{t,Z), s>0,Z{0)=t (4.8) 

(By assumptions (4.6), (4.7), it follows that t — )■ Z{t,^) is well defined on 
[0,oo).) 

Hence, for each t G [0,r], B{t) generates a Co-group {e^^^*'>)seR on L^(^) 
which is given by 

(e^WV)(0 = fiZ{s,^)), V/ G L\ff), s G R. 

It is also clear that e^^^^^V C V for all s > 0. 

Remark 4.1. Assumptions (4.5)-(4.7) can be weakened to discontinuous 
multivalued mappings ^ — > b{t,S^) satisfying (4.6), (4.7) and such that the 
solution Z = Z(s,^;t) to the characteristic system (4.8) is differentiable in 
t. The details are omitted. 

The corresponding random differential equation ()1.2p has the form 

' ||_e/3WB(t)di^^(^(^^V5g"^W^Wy) 

+ e'^B(t)^(t)e-^s(t)yds = 0, in(0,r)x^, (4.9) 

y(0,e) =a;(0° in ff, 
^ y{t,O = on {0,T)xdff. 

Then by theorem 2.1 we have 

Theorem 4.1. There exists a solution X to (j4.ip such that P-a.s. X G 
LPi(0, T; Wq'^' {&)) n LP2(0, T; VF-1'P2 (^)) n L°°(0, T; L2(^)). 
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We also note that in line with Theorem 12.21 if X„, n G N, are solutions to 
equations 

' dXl' - div5a"(t, V^X^) dt - ^bn{t, C) ■ V (6„(t, O^D dt = 
= bn{t,0 ■ V^X^d(3''{t), {t,0 in (0,r) X €?, 



Xl} = xm ff, 
XJ} = on (0, T) X 



(4.10) 



where bn ^ b uniformly on [0,T] x and an(t,y) a(t,y), a^^{t,y) 
a-\t,y), I3n{t) f3{t) for ah y € M'^, t G [0,r), then X" ^ X weakly in 
LPi(0,r; Wo'^i(^)). Standard examples refer to structural stability, PDEs 
as well as to homogenization type results for equation (|4.1|) . In latter case 
an{t,z) = a(t,nz) where a{t, •) is periodic (see e.g., [2].) 

Equation (|4.1|) is relevant in the mathematical description of nonlinear 
diffusion processes perturbed by a Brownian distribution with coefficient 
transport term 6(t,C) ■ V^X. 

The assumption pi >2 was taken here for technical reason required by the 
functional framework we work in and this excludes several relevant examples. 
For instance, the limit case pi = 1 which corresponds to the nonlinear 
diffusion function a{t, y) = Pj^, P > 0, which is relevant in material science 
and image restoring techniques (see e.g. O Sj) is beyond our approach 
and requires a specific treatment (see also [9] for the treatment of a similar 
problem with additive and continuous noise.) 

In 2-D the appropriate functional setting to treat such a problem is V = 
BV{ff) the space of functions with bounded variation on with the norm 
99(y) and H = L^i?). Here ^{y) = \\Dy\\ + J^^, |7o(y)|d^, V^V, \\Dy\\ is 
the variation of y G y, 70 (y) is the trace on di? and d^ is the Hausdorff 
measure on di?. We recall that the norm ip is just the lower semicontinuous 
closure of the norm of Sobolev space Wq'&) (see e.g., [H pag. 438].) Then 
the approach developed in section 3 can be adapted to present situation 
though V is not refiexive. We expect to treat this limit case in a forthcoming 
work. (On these lines see also [B].) 

4.2. Linear diffusion equations with nonlinear Neumann boundary 
conditions. Consider the equation 

dXt - AXt dt - ^b{t, • V^{b{t, • V^Xt) dt = 

6(t,e)-V^Xjd/3(t) in[0,r]x^ ^^.^^^ 

^Xt + C{t,Xt)BO on[0,T]x9€? 
Xo = X in ^ 

where C(*)'') = djo{t,r), \/t G (0,T), r G M and jo(t, •) is a lower semicon- 
tinuous convex function on R such that 

7i+ai|y|2 <io(t,y) <72 + a2|yP, Vy GM,t G (0,r) 
and Oj > 0, 7j G M, i = 1, 2. 

Assume also that (j4.3p holds and that b = b{t, •) satisfies conditions (j4.5p - 
(BH). 
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Then we may apply Theorems EH [221 and ESI where V = H^{ff), H 
LH^) and 



It follows so the existence of a solution X G L'^{0, T; V) n W^''^{[0, T];V') to 
(|4.1ip and also the structural stability of (|4.11|) with respect to b. Problems 
of this type arise in thermostat central. In this case 



{ai{t)H{y)+a2{t)H{-y))y if y / 
j-a2(t),ai(t)] ify = 

where Oj > 0, Vt G [0, T[ and H is the Heaviside function. 

4.3. Nonlinear stochastic porous media equation. Consider the equa- 
tion 

dXt - A^q^it, Xt) dt-\ bit, ■ V^{-Ay\b{t, • V5((-A)"iXt) dt 
= 6(t,0 • V5(-A)-iXtd/3(t), in (0,r) X ^ 

Xq = X in ^ 
Xt = on (0, oo) X dff 

(4.12) 

Here 0' C R'^, d = 1,2,3 is a bounded open domain and (—A) ^ is the 
inverse of the operator Aq = -A, D{Ao) = H^{^) D H'^{ff). The function 
(f): (0,T) X R'^ — 7- M is assumed to satisfy the following conditions 

(k) (j) = (f){t, r) is monotonically decreasing in r, measurable in t and its 
potential 

j{t,r) = / (t>{t,T)dT, te{0,T) 
Jo 

satisfies the growth conditions 
71 + ai |r|Pi <j(t,r) < 72 + a2|r|P2, Vr G M, a; G f^, t G [0, T] (4.13) 

j{t,-r) <cij{t,r)+C2, VrGM,tG(0,r) (4.14) 
where ^ < pi < p2 < oo if d = 3, 1 < pi < p2 < oo if d = 1,2. 

Then equation (|4.12p can be written as (jl.ip . where H = H^^{ff), V = 
LP^{ff) and A{t) = d^{t, ■) where ■): H ^Ris defined by 

/ j{t,y)di ifyeH-\i?),j{t,y)eLU 

+00 otherwise, 

and B{t),teW^ is defined by 

B{t)u = b{t, i) ■ V((-A)-^'u), u£V. (4.15) 

The space V' is in this case the dual oi V = LP^ii^) with H~^{ff) as pivot 
space. By the Sobolev embedding theorem it is easily seen that since pi > | 
we have V C H~^{ff). The scalar product on H is defined by 

{u,v)H-i(ff) = u{z), z = {-A)~^v. 

It is well known that A(t)X = —A^(f){t, X) is indeed the subdifferential of 
ip{t,-) in i7"i(^) (see e.g., H pag. 68]). 
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As regards b: [0,T] x ^ R"^ we assume that conditions (|i3|) ~- (|i7ri) hold. 
We note that for each t G [0,r], B{t) E L{V,H-^{ff)) is densely defined 
and skew-symmetric on H^^{&) = H. Indeed we have 

{B{t)u,u)= [ div(6(t,e)(-A)-iu) • (-A)-^ud^ = 
Jff 

= \ [ 6(t,o-v|(-A)-MOI'rf^ = o, 

because div^fe = and b{t,^) ■ = on dff. 

Moreover, for each t € [0,T], B[t) is rn-dissipative on H~^{ff). Indeed for 
each / G H~^{ff), the equation u — B{t)u = f can be equivalently written 
as V = {—A)~^u, where 

-Av - b{t, ■) -Vv = f in ff, 
V = on Off. 

By Lax-Milgram lemma the latter has a unique solution v G H^^ff) and 
therefore u G H~^{ff) as claimed. Moreover, if / G L^^ (ff) and dff is of class 
then by the Agmon-Douglis-Nirenberg theorem v G W'^'^^ {ff)r\WQ^\ff) 
and so u eV. 

Hence, B{t) generates a Co-group {e^^^^^)s^R on H = H^^{ff) which leaves 
V = LP^{ff) invariant. 

Then we may apply Theorem 12.11 as well as the approximation Theorem 12.21 
to the present situation. We obtain 

Theorem 4.2. There is a unique solution X to (j4.12|) such that P-a.s. X G 
LP^{0,T;LP^{ff))nL'^{0,T;H'^{ff)). Moreover, the solution X is a limit 
of approximating solutions when the Brownian motion (5 is approximated by 
a sequence of smooth processes. 

Moreover, \l (pn ^ <P and (p^ —?■ (p* , bn ^ b we find by Theorem 2.2 that 
the corresponding solutions Xn to ()4.12p are convergent to solution X to 
(|4.11|) . The details are omitted. 

Existence for stochastic porous media equation of the form 

dXt - A^^{Xt) dt = a{Xt) dW{t) in (0, T) x ff 
Xo = x in ff 
Xt = on (0, T) X dff, 

when Wt is a Wiener process of the form 

oo 

W{t,0 = Y.l'kek{OMt) 

k=l 

with X^fcLi Aifc A| < oo, Acfc = — A^Cfc in ff, e^ G HQ{ff), and a = a{x) 
is a linear continuous operator, were studied in [HI [9]. Note that in this 
case the noise term can also be written in our form with commuting and, 
contrary to our paper, bounded operators Bj. Here the multiplicative term 
a{Xt) = b ■ V{—A)~^Xt is however discontinuous on the space H~^{ff) and 
so Theorem 4.2 is from this point of view different and in this sense more 
general. 

Equation (j4.12p models diffusion processes and the motion of fluid flows 
in porous media. The case considered here {pi > 1) is that of slow diffusion. 
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Remark 4.2. Theorem 12.41 and Remark 12.11 are also valid in the current 
setup. 

Appendix A. Convex functions 

We summarize in this paragraph some facts about convex functions, which 
we have used in our paper. 

Given a convex and lower-semi continuous function : y — )• M = (— oo, oo] 
we denote hy d(j) : Y ^ Y' (the dual space) the subdifferential of (j), i.e. 

dcPiy) := {zeY': ^ - 0(n) < {y - u, z) , eY] . (A.l) 

(Here (•, •) is the duality paring between Y and Y'). The function (f)* -.Y' ^ 
Y defined by 

<l)*{z) = sup{(2/,z)-0(y) :yey}, (A.2) 

is called the conjugate of (f) and as (j). Similarly to it is convex lower semi con- 
tinuous function on Y' . Also we notice the following key conjugacy formulae 
(see e.g. 13 p. 89]). If ?/ G Y, and z G 

m + (t^*{z) > {y,z) (A.3) 

^{y)+^*{z) = {y,z) iSzedcbiy), (A.4) 

A vector x* is said to be a subgradient of a convex function cj) at a point 
X if 

(t>{z) > (t>{x) + {x* , z - x) . (A. 5) 

Moreover, straightforward calculations give 

^*{x*) = cPl{x*)-{y,x*), (A.6) 
whenever (/)(x) = (l)y{x + y). 
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